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■ ABSTRACT 

f-n | The Planck Australia Telescope Compact Array (Planck- ATC A) Co-eval Observa- 

tions (PACO) have provided multi-frequency (5-40 GHz) flux density measurements 
I 1 of complete samples of Australia Telescope 20 GHz (AT20G) radio sources at fre- 

quencies below and overlapping with Planck frequency bands, almost simultaneously 
with Planck observations. In this work we analyse the data in total intensity for the 
spectrally-selected (SS) PACO sample, a complete sample of 69 sources brighter than 
<S20GHz = 200 mJy selected from the AT20G survey catalogue to be inverted or up- 
turning between 5 and 20 GHz. We study the spectral behaviour and variability of 
the sample. We use the variability between AT20G (2004-2007) and PACO (2009- 
2010) epochs to discriminate between candidate High Frequency Peakers (HFPs) 
and candidate blazars. The HFPs picked up by our selection criteria have spectral 
peaks > 10 GHz in the observer frame and turn out to be rare (< 0.5% of the 
•S^oghz 200 mJy sources), consistent with the short duration of this phase implied 
by the 'youth' scenario. Most (~ 89%) of blazar candidates have remarkably smooth 
spectra, well described by a double power-law, suggesting that the emission in the 
PACO frequency range is dominated by a single emitting region. Sources with peaked 
PACO spectra show a decrease of the peak frequency with time at a mean rate of 
— 3±2GHzyr~ 1 on an average timescale of (t) = 2.1±0.5yr (median: r mc dian = l-3yr). 
The 5-20 GHz spectral indices show a systematic decrease from AT20G to PACO ob- 



servations. At higher frequencies spectral indices steepen: the median afg is steeper 
than the median by da = 0.6. Taking further into account the Wide-field Infrared 
1 Survey Explorer (WISE) data we find that the Spectral Energy Distributions (SEDs), 

vS(v), of most of our blazars peak at ^p ED < 10 5 GHz; the median peak wavelength 
is Ap ED ~ 93 fim. Only 6 have j/| ed > 10 5 GHz. 

Key words: galaxies: active - radio continuum: galaxies - radio continuum: general 
- cosmic microwave background. 



1 INTRODUCTION 

The PZancfc-Australia Telescope Compact Array (Planck- 
ATCA) Co-eval Observations (PACO) project consists in 
observations with the Australia Telescope Compact Array 
(ATCA) of three complete samples drawn from the Aus- 
tralia Telescope 20 GHz (AT20G) catalogue (Murphy et 
al. 2010; Hancock et al. 2011) in three pairs of 2GHz-wide 
bands (5.5 and 9 GHz, 18 and 24 GHz, 33 and 39 GHz). The 
selected samples are: the 'bright PACO sample' including 
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sources with 520GHz ^ 500 mJy (Massardi et al. 2011b); the 
'faint PACO sample' (Bonavera et al. 2011) going down to 
S20GHZ = 200 mJy in the South Ecliptic Pole region; and the 
'spectrally-selected PACO sample', presented in this paper. 
The PACO observations were carried out in several epochs 
between July 2009 and August 2010. At least one obser- 
vation was made within 10 days from the Planck satellite 
observations in any of the LFI channels (30, 44, 70 GHz). 
Observations and data reduction procedures have been pre- 
sented in Massardi et al. (2011b). The comparison between 
PACO data for all the samples and measurements in the 
Planck Early Release Compact Source catalogue (ERCSC; 
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Planck Collaboration VII 2011) will be presented in Mas- 
sardi et al. (in preparation). 

A particularly interesting class of sources with spec- 
tra rising at frequencies above 5 GHz are the 'High Fre- 
quency Peakers' (HFPs; Dallacasa et al. 2000), believed to 
be the most compact and youngest members of the Giga- 
hertz Peaked Spectrum (GPS) source family (e.g., O'Dea 
1998; Tinti & De Zotti 2006). The first HFP selections (Dal- 
lacasa et al. 2000; Stanghellini, Dallacasa, & Orienti 2009) 
had peak frequencies in the observer frame mostly in the 
range 5-10 GHz. The AT20G survey offered the possibil- 
ity of selecting complete samples of more extreme HFPs, 
with spectral peaks above 10 GHz (Hancock et al. 2010). 
The well established correlation between intrinsic peak fre- 
quency, Vm, and projected linear size (O'Dea 1998) indicates 
that v m ~ 10-20 GHz corresponds to scales ~ 10 pc. Multi- 
epoch VLBA and VLBI observations of Compact Symmetric 
Objects (CSOs), as genuinely young HFP/GPS sources are 
generally found to be (Orienti et al. 2006; Orienti & Dalla- 
casa 2012), showed sub-relativistic expansion velocities, fre- 
quently of 0.1-0. 2c (Gugliucci et al. 2005; Polatidis 2009; 
An et al. 2012). Thus the AT20G HFPs are expected to 
have ages of a few hundred years and therefore to give us an 
insight on the initial stages of radio activity. The available 
data indicate that these early phases are characterised by 
high luminosity, fast evolution in size and peak frequency 
(Begelman 1996; Dallacasa 2003; Tinti & De Zotti 2006), 
and must therefore be very rare. Only very large area, high 
frequency surveys, such as the AT20G that covers the en- 
tire Southern sky, can thus provide statistically significant 
samples or constraints. 

A peaked spectral shape, however, does not guarantee 
that a source is truly young. Beamed radio sources (blazars) 
can also show peaked spectra when a flaring synchrotron 
component, self-absorbed up to high frequencies, dominates 
the emission spectrum. But the latter sources change their 
spectral shape keeping convex spectra for a relatively short 
time and may show a flat spectrum when re-observed af- 
ter a few years, and anyway exhibit substantial variability 
(Tinti et al. 2005; Torniainen et al. 2007). In contrast, truly 
young sources maintain their convex spectra and show low 
or undetectable variability, apart from the spectral evolu- 
tion associated to the source expansion (Orienti, Dallacasa, 
& Stanghellini 2010). Thus multi-epoch, multi-frequency ob- 
servations of peaked spectrum sources proved to be very ef- 
fective in discriminating among the two populations. We use 
the variability properties between AT20G (2004-2007) and 
PACO (2009-2010) epochs to discriminate between young 
sources and blazars. 

Our search for HFP candidates differs from the ear- 
lier study by Hancock et al. (2010), also based on the 
AT20G survey, in several respects. We start from a com- 
plete flux limited sample selected over the whole South- 
ern sky, while the sample by Hancock et al. has no well 
defined flux density limit and is restricted to the region 
18 h < RA < 01 h 30 m . Also, their variability threshold for 
HFP candidates is mostly based on data at a single fre- 
quency (20 GHz) on relatively short timescales (1-3 yr), 
while we have multi-frequency data on somewhat longer 
timescales (typically 2-4 yr, but up to 6 yr). On the other 
hand, Hancock et al. obtained radio observations also at 95 
GHz (while PACO observations span the range 5-40 GHz) 



as well as optical imaging and spectroscopy of a subset of 
their sources. As for blazars, the combination of AT20G and 
PACO data can constrain the duration of high frequency 
flares and the distribution of their peak frequencies. 

The structure of this paper is as follows. In §[2] we de- 
scribe the data for the spectrally-selected PACO sample. 
In §[3] and in §[4] we discuss candidate HFPs and candidate 
blazars, respectively. Finally, in §[5] we summarise our main 
results. 



2 SAMPLE CHARACTERISATION 
2.1 Sample selection 

The spectrally-selected (SS) PACO sample comprises all the 
69 sources with &oghz > 200 mjy and spectra classified in 
the AT20G catalogue (Massardi et al. 2011a) as inverted 
or upturning in the frequency range 4.8-20 GHz, selected 
over the whole Southern sky. These sources have been re- 
observed for the PACO project between September 2009 and 
February 2010, with a scheduling process optimised to allow 
observations at all the frequencies almost simultaneous (i.e. 
within 10 days) with the Planck satellite. Fifteen SS sources 
are also part of the faint PACO sample (5*20GHz > 200 mjy, 
Bonavera et al. 2011) and 21 of the bright PACO sam- 
ple (S20GHZ > 500 mjy, Massardi et al. 2011b), including 
6 sources in common with the faint sample. 

The ATCA observations of each source over the fre- 
quency range 5-40 GHz have been carried out in the 
same day but, due to bad weather conditions, in one case 
(J173340- 793555) a complete spectrum for the source has 
been obtained by joining the observations done in two sub- 
sequent days. The spectra of two other sources (J040019- 
225624 and J184827- 735337) are not reliable as the sources 
turned out to be extended while the flux estimation relies 
on the point-like assumption. Data on these two sources are 
presented in the catalogue but not included in the analysis, 
so hereafter we consider a sample of 67 sources. 

The classification of Massardi et al. (2011a) has been 
done by considering two spectral indicef]: a|, between 4.86 
and 8.56 GHz, and , between 8.56 and 19.9 GHz. At the 
time of the AT20G observations the sample was defined by 
having a| > and af > ("inverted") or «f < and 
al° > ("upturning"), after having excluded the sources 
with both < 0.5 and |aj°| < 0.5 ("fiat"). 

We have redone this classification using the PACO mea- 
surements over the 512 MHz-wide sub-bands centred at the 
frequencies nearest to the AT20G ones (4.732, 8.744 and 
18.768 GHz). The best quality PACO measurements, i.e. 
those referring to the epoch at which the highest accuracy 
was achieved (due to the combination of good weather, mini- 
mal radio interferences and instrumental malfunctions, etc.), 
were selected for this purpose. As illustrated by Fig.fT]a large 
fraction of the sample has significantly modified its spectral 
behaviour, moving to the "flat-spectrum" region of the di- 
agram. This already shows that the majority of sources in 
the SS sample are blazars, caught in a bright stage by the 
AT20G survey, that decayed to a more quiescent phase at 
the time of PACO observations, i.e. on a timescale of a few 

1 We adopt the convention S„ <x u a 
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Figure 1. Colour-colour diagram (a| versus ct|) of the sources 
in the SS PACO sample. Green triangles: AT20G data on which 
the sample selection was made; black diamonds: PACO data. The 
dotted lines show the q| = and Og = axes and the region 
defined by |a|| < 0.5 and < 0.5 contains sources classified 

as 'flat-spectrum' by Massardi et al. (2011a). 



to several years. We will come back to this in Sec.[3l when 
we will discuss the selection of "genuine" HFP sources. 



2.2 Spectral behaviour 

As a first step we have checked whether the spectra of our 
sources can be adequately described, over the full frequency 
range 5-40 GHz, by the double power-law 



S(y) = So/[(u/va)-° + {v/v y b ], 



(1) 



used by Massardi et al. (2011a) and Bonavera et al. (2011) 
for the bright and faint PACO sample, respectively. So, vq, 
a and b are free parameters. Such a smooth spectrum may 
not necessarily fit the data for blazars whose radio emission 
results from the superposition of peaked spectra of several 
parts of the jet racing towards the observer at highly rel- 
ativistic speeds and self-absorbed up to different frequen- 
cies (e.g. De Zotti et al. 2010). The fit has been performed 
in logarithmic units by minimising the x 2 with a nonlinear 
optimisation technique based on an implementation of the 
Generalised Reduced Gradient optimisation method. Unless 
otherwise noted, the results reported below refer to the "best 
quality" PACO spectra (see 5 l2~Tj) . 

We have considered as bad fits those with \ 2 ~ (x 2 ) > 
3u x where (x 2 ) and a x are mean and standard deviation 
of the Gaussian fit of the reduced x 2 distribution over the 
whole PACO sample. Those are respectively (x 2 ) = 1-3 and 
a x — 0.54. We note that, in general, the x 2 values associated 
to the best fits are relatively high. This is due to the fact 
that, over the wide frequency range considered, the observed 
spectra show some small roughness, detected by our high 
accuracy measurements, so that a smooth function cannot 
provide a perfect fit. Still, the double power-law fits look 
acceptable for most sources (see Figs. [2] [6]). 



Frequency 


PACO-AT20G 


PACO-PACO 


[GHz] 


(2-4 yrs) 


(6 months) 


5 


13.1 ± 1.9 


5.5 ± 1.5 


9 


15.0 ± 2.2 


4.9 ± 1.7 


18 


20.3 ± 2.5 


5.7 ± 1.7 


24 




8.2 ±2.0 


33 




4.8 ±3.0 


39 




5.3 ±3.2 



Table 1. Median variability indices [%] of the SS sources and 
associated errors (1.25 cr/^/N — 1, where N and a are the num- 
ber of the variability indices and the standard deviation of their 
distribution, respectively). 



For the SS sample, the double power-law model has 
been accepted for 57 out of 67 sources, i.e. in 85% of the 
cases. For them we have computed low- (5-10 GHz) and 
high-frequency (30-40 GHz) spectral indices based on the fit 
and found that they fall in the following spectral categories: 

• "peaked": af > 0.3, af < -0.3; 20 sources (30%) 

• "down-turning": afpj ^ min(a5°— 0.35, — ajj°) and al° ^ 
0.3; 26 sources (39%) 

• "self-absorbed": al° 
7 sources (10%) 

• "flat" with la! 1, |a|g| «S 0.3: 4 sources (6%). 

No source switched from "inverted" or "upturning" accord 



- 0.35 ^ af , > max(-0.3, -af 



ing to the AT20G data to "steep" (a 5 u , < -0.3) ac- 
cording to the PACO observations. Nor any source kept 
an "inverted" spectrum (Q5 , a|o > 0.3) up to the high- 
est frequencies in the PACO range. A similar result was 
found for the bright and faint PACO samples. The peak 
frequency v p can be derived from the fitting parameters as 
^ = ^o(-6/a) (1/(6 " a)) - 

The remaining 10 (15%) sources have complex spectra 
(see §U): besides a main peak, they also exhibit a sort of 
plateau, which can either be at v lower (6 sources) or higher 
(4 sources) than the peak frequency. This behaviour could 
be indicative of multiple emission components. 



2.3 Variability 

Following Sadler et al. (2006) we define the variability index 
at the frequency v as: 



V 



100 k /E[^-(^)] 2 -E«) 



{S") 



(2) 



where Si and 0" are the flux density and the associated 
error of a given source measured at the i-th epoch, n is the 
number of epochs, and (S v ) is the mean of the available 
measurements. For 10-15% of the sources (depending on 
frequency and time lag) y^fS*" — (S v )] 2 < E(o"j") 2 , so that 
the argument of the square root in eq. (f2I is negative. In 
these cases we estimated the 1 a upper limit to the variability 
index of the source as: 



100 /EK- 



(3) 



The comparison between AT20G (2004-2008) and 
PACO (2009-2010) measurements gives information on the 
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Figure 2. PACO data (crosses) for the sources identified as HFP, with different colours corresponding to different epochs. The dotted 
lines are the spectral fits the to the PACO data for the epoch with the most accurate data. Also shown for comparison are the AT20G 
data (triangles). 



source variability at ~ 5, 8, and 20 GHz (Kms> Kms> Kms) 
on typical timescales of ~ 2-4 years (but of up to 6 years in 
some cases), while the comparison of PACO measurements 
at different epochs informs on variability on a timescale of 
~ 6 months, since they were coupled to the Planck scans 
which have a ~ 6 month periodicity. The variability indices 
on the latter timescale have been computed, for all PACO 
frequencies, selecting for each source the best (smallest error 
bars) pair of observations spaced by 6 months within ±20%. 

The median PACO-AT20G and PACO-PACO variabil- 
ity indices and the associated errors are reported in Table[T] 
The PACO-AT20G variability of the SS sample turns out 
to be systematically higher than that found for the bright 
PACO sample (11.8 ± 2.1, 11.7 ± 2.2 and 13.9 ± 1.9 for 5, 
8 and 20 GHz respectively; Massardi et al. 2011b), as ex- 
pected since the SS sample contains a larger fraction of flar- 
ing objects. There is also a trend towards higher variability 
at higher frequencies, also found for the bright PACO sam- 
ple (Massardi et al. 2011b). On the other hand the PACO- 
PACO variability amplitude is similar to that found for the 
bright sample on the same timescale (6 months). No trend 
with frequency shows up in this case. Of the three spec- 
tral subclasses containing at least 10 objects, the "peaked" 
one turns out to be the least variable. On the 6 months 
timescale their median variability index is in the range 1.6- 
7.6, depending on frequency; the median variability index of 
"down-turning" sources ranges from 6.4 to 15.0 and that of 
"complex spectrum" sources ranges from 4.7 to 17.9. On the 
AT20G-PACO timescale (2-4 yr) the median variability in- 
dices are in the range 9-10, 13-32 and 15-40 for "peaked", 
"down-turning" and "complex spectrum" sources, respec- 
tively. 

For 15 out of the 20 "peaked" sources we could de- 
termine the timescale for the variation of v v in the source 
frame, r = v p / (Av p / At Icst ). Here v p is the peak frequency 
at the first epoch at which it could be determined, Av p is 
its variation after a time Ai b s , in the observer frame, and 
At res t = At ba/(1 + z) is the corresponding time lag in the 
source frame. The redshift z was available for 3 out of the 
15 sources; for the remaining ones we used the median value 
over the candidate blazars of the sample having a redshift 
estimation (z — 1.33). There is a trend towards a decrease of 
i/ p with time: Av p /At les t is negative for 12 of the 15 sources; 



its mean value is — 3 ± 2GHzyr _1 and its median value is 
— 2.9GHzyr _1 . This trend is expected if the size of the flar- 
ing region expands, and correspondingly its self-absorption 
frequency decreases, as it flows along the jet. The average 
variation timescale is (r) = 2.1 ± 0.5 yr and the median is 

^"median — 

1.3 yr. There is no detectable correlation between 
r and v v \ this is not surprising since a correlation would 
anyway be hard to detect given the limited range of peak 
frequencies of our sources. 



3 SELECTION OF HFP CANDIDATES 

As mentioned in §[T] variability studies are a powerful tool 
to distinguish GPS (and HFP) sources from flaring blazars. 
GPS sources should typically vary by no more than 10% 
(O'Dea 1998) while blazars are in general much more vari- 
able. Six objects (out of the 67 included in the sample, i.e. ~ 
9%) are found to have V" ms ^ 10%. However, a low variabil- 
ity on the limited timescales of our observations is not a suf- 
ficient condition for sources not to be blazars, since blazars 
may stay in a quiescent state for several years. In fact, 
additional data taken from the NASA/IPAC Extragalactic 
Database (NED) and from the Massaro et al. (2009) blazar 
catalogue do show evidence of strong variability for 3 out of 
the 6 sources, that therefore are classified as blazars. Thus 
only 3 HFP candidates remain (J045316-394905, J081240- 
475153, J183923-345348), i.e. 4.5% of sources in the sample. 
They are all close to the 20 GHz flux density threshold. We 
list their relevant properties in Table [2] while their spectra 
are shown in Fig. [2] Their peak frequencies, v p , in the ob- 
server's frame, are in the range 15-22 GHz. 

There is a well established correlation between the rest 
frame peak frequencies, v m = v v (\ + z), and the projected 
linear size I (O'Dea & Baum 1997) 

log(i/ m /GHz) ~ -0.21(±0.05) - 0.65(±0.05) log(Z/kpc), (4) 

implying that intrinsic peak frequencies > 10 GHz corre- 
spond to scales of a few tens of parsecs or less. As mentioned 
in §[T] GPS sources are generally found to have expansion 
velocities of 0.1-0. 2c so that sizes of ~ 10 pc correspond to 
ages of hundreds of years. 

It must be noted that the selection criteria for 
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Table 2. HFP candidates. 



AT20G name 


RA [deg] 


DEC [deg] 


«5° 


Q 30 


u p [GHz] 


z« 


v m [GHz] 


J[pc] 


J045316-394905 


73.3191 


-39.8182 


0.61 


-0.16 


15.66 


0.6 


24.3 


3.5 


J081240-475153 


123.171 


-47.8647 


0.50 


-0.45 


15.45 


0.7 


26.4 


3.1 


J183923-345348 


279.848 


-34.8968 


1.39 


-1.05 


21.62 


0.5 


32.2 


2.3 



M Estimated from B magnitude from superCOSMOS as in Burgess & Hunstead (2006). 
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Figure 3. PACO data (crosses) for the sources classified as spectrally peaked blazars with different colour corresponding to different 
epochs. The dotted lines are the spectral fits to PACO data for the 'best' epoch. Also shown for comparison are the AT20G data 
(triangles). 



GPS/HFP sources are somewhat fuzzy. For example, the 
sample of 21 GPS/HFP candidates selected by Hancock et 
al. (2010) includes 5 sources with S20GHZ > 200 mjy not in- 
cluded in our sample because they have both a§ < 0.5 and 
a|° < 0.5, and were therefore classified as "flat-spectrum" 
by Massardi et al. (2011a). Three out of these 5 candidates 
were found to have higher variability than expected for gen- 
uine GPS/HFP galaxies, one was classified as a likely and 
one as a possible GPS/HFP galaxy. This suggests that our 
HFP sample may be incomplete. Since the area covered by 
the Hancock et al. sample is about 30% of the AT20G area, 
the incompleteness may amount to a factor of ~ 2. On the 
other hand, the amount of variability for the Hancock et al. 
sample has been mostly determined on the basis of observa- 
tions at only two epochs with time lags of 1-3 years and only 
at 20 GHz. As pointed out by Hancock et al. themselves, well 



monitored AGNs show only small flux density variations on 
1-2 year timescales, with larger outbursts occurring roughly 
once every 6 years. Hence more data are necessary to assess 
the GPS/HFP nature of the two additional candidates. 

On the whole, we may conclude that the fraction of 
HFPs in our complete sample of inverted-spectrum sources 
with 5*20GHz > 200 mjy, albeit uncertain, is < 10%. Thus 
the high-frequency peaking population is largely dominated 
by relativistically beamed objects (blazars). A similar con- 
clusion was reached by Bolton et al. (2006) for a sample 
of 9C sources brighter than 25 mjy at 15 GHz with spec- 
tra peaking at u v > 5 GHz. Since SS sources are only 4.9% 
of the 1360 AT20G sources brighter than that flux density 
limit, the fraction of HFP candidates among AT20G sources 
is < 0.5%. The rarety of genuine HFPs is indicative of a very 
short duration of the HFP phase, as expected in the frame- 
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Table 3. Blazar candidates in our sample. The blazar type is taken from Massaro et al. (2009) for the 24 sources in BZCAT. The spectral 
type is based on PACO data (see § 12.21 1. 



AT20G name 


RA [dog] 


DEC [deg] 


blazar type 


redshift 


spectral type 


up [GHz] 


v m [GHz] 


J UUZOlO-OOlZ'rty 


O.OOoZo 


-00.Z100 






peaked 


iy .00 
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Figure 4. Same as Fig. [3] for the sources classified as "down-turning" blazars. 



work of the 'youth' scenario for these sources. For example, i.e. of hundreds of years. According to the model by Snellen 
the model by Begelman (1996) implies u m oc r~ Se (Tinti & et al. (2000) u m oc r ~ 5 / 18 an d the evolution timescale of v m 
De Zotti 2006) where r is the source age, 5 ~ 0.65 is the is only moderately larger, 
slope of the v m -svz.e relation [eq. @] and e — 3/(5 — n), n 
being the slope of the radial density profile (p e oc r~ n ) of the 
medium surrounding the expanding source. For the values of 
n (n = 1.5-1.9) favoured by Begelman (1996), the evolution 
timescale of v m is within a factor of two of the source age, 
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Figure 5. Same as Fig.[3]for the sources classified as "self-absorbed" blazars. 
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Figure 6. Same as Fig. [3] for the sources classified as "flat" blazars. 
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Figure 7. Same as Fig[3]for the sources identified as blazars and classified as "complex spectrum". A polynomial fit is shown instead of 
the double power-law one (dotted line). 
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Figure 8. Left: distribution of (hatched histogram) and a|JJ (filled histogram) spectral indices for the sources with double power- 
law spectra, compared with the distribution of AT20G a|° for the same sample of sources (empty histogram). Right: distribution of 
peak frequencies for SS sources with double power-law spectra (peaked, down-turning and self-absorbed) in the observer's frame (filled 
histogram) and in the rest frame for the subset of sources having a redshift measurement or estimate (hatched histogram). 
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Figure 10. Left-hand panel: distribution of the peak f Pi sED obtained by fitting a parabolic model to the PACO data combined with 
WISE data. The > signs correspond to the 6 objects for which only a lower limit to fp BD could be obtained. The < signs correspond 

to objects missing WISE photometry, for which we could obtain only upper limits to ^ ED . Right-hand panel: peak luminosity versus 
,.SED 



4 BLAZARS 

All SS sources not classified as HFP candidates (64 sources, 
95%) are likely blazars. A positional cross-match with the 
Third Edition (April 2011) of the multi-frequency blazar 
catalogue (BZCAT) by Massaro et al. (2009) allowed us to 
confirm the blazar classification for 24 of them. The list of 
sources with coordinates and spectral properties is given in 
Table O The spectral classification and the value of v v are 
based on PACO data as described in Sec. 12.21 The blazar 
type is from the BZCAT where blazars are subdivided into 
flat-spectrum radio loud QSOs (QSO RLoud), BL Lacs and 
Uncertain type. With only one exception (Uncertain type) 
all our blazars included in BZCAT (which are however a 
minority, ~ 40% of the total number of our blazars) are 
classified as QSO RLoud. 

The spectra based on PACO observations are shown in 



Figs. |3HZ1 (one figure per spectral type) for all epochs, to- 
gether with AT20G (triangles) data points. Even though 
the sources were selected for having inverted or upturn- 
ing AT20G spectra, the PACO observations show, for most 
of them, fiattish spectra in the 5-20 GHz range. The me- 
dian PACO spectral index a|° is 0.14, to be compared with 
the median AT20G value of 0.48 for the same sources (see 
the left panel of Fig. [8}. This indicates that these sources 
were caught by AT20G observations during a bright emission 
phase of a dense region probably in the inner part of the jet, 
and that the emission decayed by the time of PACO follow- 
up. The spectra show a significant steepening at higher fre- 
quencies (see again the left panel of Fig. [5} : the difference 
between the median a|° and the median 

alo is 0.6. Only 

one source, J040353-360500, has af Q > (afg = 0.26); its 
spectrum however flattens towards the top of the frequency 
range (afg = -0.006). 
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Figure 9. Representative examples of SEDs including PACO and 
WISE data. Upper row, from left to right: two "down-turning" 
and a "self-absorbed" blazar. Lower row, from left to right: a 
"peaked", a "flat" and a "complex spectrum" blazar. 



The distribution of peak frequencies of SS sources with 
double power-law spectra is displayed in the right-hand 
panel of Fig. [8] The mean peak frequency of blazars, in 
the observer's frame, is 14 GHz with a standard deviation 
of 6 GHz. The candidate HFPs lie in the low- frequency part 
of the distribution. The median rest-frame peak frequency 
of the 24 blazars with redshift information is 20 GHz. 

Interestingly, with only one exception, the 'best epoch' 
spectra of sources for which a double power-law [eq. (0] 
fit could be obtained have v v < 40 GHz. The exception is 
J040353-360500 (see Fig. © for which Planck ERCSC data 
are available; a double power-law fit including them yields 
v p ~ 51 GHz. 

Since all SS sources have high-frequency spectral in- 
dices ail > -1 their SEDs vS(y) are rising in the frequency 
range covered by PACO observations. To determine the SED 
peak frequency, fp ET> , we constrained the descending part of 
the spectrum with WISE (Wright et al. 2010) data, at 3.4, 
4.6, 12 and 22 tun (~88000-13600 GHz). Of the 64 blazars, 
51 (79%) have a counterpart in WISE; for the remaining 13 
(21%) we used the upper limits given by Wright et al. (2010): 
0.08,0.11, 1 and 6 mjy at 3.4, 4.6, 12 and 22 (jm respectively. 
We fitted the PACO plus WISE data with a parabola in the 
log i/-log uS(u) plane, parametrised by the th , 1 st and 2 nd 
order coefficients Co, cr and C2, and got Vp EU = — Ci/(2c2). 
The accuracy of this estimation is limited by the gap be- 
tween the PACO and the WISE frequency ranges and by 
variability, as PACO and WISE data are not simultaneous. 
In two cases the PACO plus WISE data turned out to be 
incompatible with a parabolic fit, but for the other objects 
the fits looks reasonably good. Some representative exam- 
ples are shown in Fig. [5] 

The distribution of Vp BU is shown in the left-hand panel 
of Fig. [10] Most of them (56) have ^ ED < 10 5 GHz and 
are therefore classified as "Low Synchrotron Peak" (LSP) 
according to Abdo et al. (2010). Their median peak fre- 
quency (log(^ p /GHz) = 3.51) corresponds to a peak wave- 
length Ap ED ~ 93 ^m; there is no significant trend of the 
peak frequency log(i/ p /GHz) with the PACO spectral class. 

For the remaining 6 sources (J042127-481805, J093841- 



655611, J182927-570534, J194131-760548, J203958- 720655, 
J230303-184126) the fit yields values of i/| ED larger than 
10 s GHz, i.e. above the range covered by the data. Thus, 
only a lower limit (fp ED 10 s GHz) could be obtained. 
The Two Micron All Sky Survey (2MASS; Skrutskie et 
al. 2006) photometry, available for all these sources, con- 
firms that their SEDs keep flat up to at least the J-band 
(2.4 x 10 5 GHz). B-band photometry, given in the NED 
for 3 of them (J042127-481805, J093841-655611, J230303- 
184126), is consistent with SEDs being still flat. Thus, these 
are likely "High Synchrotron Peak" (HSP) sources (^p ED > 
10 6 GHz). The sharp drop of the J/p ED distribution below 

,SED 



]3| 

SED 



Up = 10 GHz is a selection effect: sources with lower 



values of Vp are too faint to be detected by WISE and 
therefore we could only obtain upper limits ( Si 10 3 GHz) to 
their peak frequencies. 

The addition of WISE data gives us a substantially dif- 
ferent perspective on source spectra. Three of the 6 likely 
HSPs were classified as "down-turning" on the basis of 
PACO data, 2 other likely HSPs were classified as "peaked" , 
and the remaining one was classified as "complex spectrum" . 
This illustrates once more the complexity and the time- 
dependence of blazar spectra. 

The right-hand panel of Fig. \W\ shows z^p ED versus the 
peak luminosity [^p ED L(i>J! BD )]. According to the blazar se- 
quence model (Fossati et al. 1998; Ghisellini et al. 1998), 
the synchrotron peak frequency is anti-correlated with the 
radio luminosity. Our data do not show any statistically sig- 
nificant correlation: the Spearman correlation coefficient is 
p = —0.05, implying that the null hypothesis (no correla- 
tion) is rejected only at the 0.4 a level. 



5 CONCLUSIONS 

We have presented an analysis of the "spectrally-selected" 
(SS) PACO sample comprising all the AT20G sources with 
S20GHZ > 200 mjy and spectra classified by Massardi et 
al. (2011a) as inverted or upturning in the frequency range 
5-20 GHz. The sample consists of 69 sources, two of which 
were however rejected because they turned out to be ex- 
tended and thus to have unreliable flux density estimates. 
The PACO observations were carried out in three pairs of 
2 GHz-wide bands (5.5 and 9 GHz, 18 and 24 GHz, 33 and 
39 GHz) in several epochs between July 2009 and August 
2010. At least one observation was made within 10 days 
from the Planck satellite measurements. 

Most (85%) of the source spectra are remarkably 
smooth and well described by a double power-law. Twenty 
have a peaked spectrum. Most (12 out of 15) of those for 
which PACO observations allowed an accurate determina- 
tion of the spectrum at different epochs show a decrease 
of the peak frequency with time. The mean decrease rate 
for the 15 peaked sources is —3 ± 2GHzyr- 1 . The average 
variation timescale is (r) = 2.1 ± 0.5 yr and the median is 

Tmcdian = 1-3 yr. 

Only three of our sources are found to be good candi- 
date High Frequency Peakers, with peak frequencies, in the 
observer's frame, v v > 10 GHz. Their 20 GHz flux densities 
are close to the sample threshold and their peak frequencies 
lie in the lower part of the distribution for the full sample. 
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Even allowing for a generous correction for the possible in- 
completeness of our selection, these extreme (y v > 10 GHz) 
HFPs are only a small fraction (< 0.5%) of sources with 
5*20GHz > 200 mjy, consistent with the very short duration 
of the HFP phase expected in the framework of the 'youth' 
model. 

All the other sources in our sample are likely blazars, 
and 24 of them, all but one classified as flat-spectrum radio 
QSOs, are listed in the last edition of the blazar catalogue 
by Massaro et al. (2009; BZCAT). We find a clear decrease 
of 5-20 GHz spectral indices from AT20G to PACO obser- 
vations, implying that these sources were generally caught 
by the AT20G survey during a bright phase at 20 GHz. At 
higher frequencies spectral indices steepen: the median afo 
is steeper than the median a§° by 8a = 0.6. All but one of 
the sources for which a peak of the flux density distribution 
could be defined have v v < 40 . This is at odds with the 
results by Hancock et al. (2010) who found that 12 of their 
21 sources have v v > 40 GHz, including 4 sources peaking 
above 95 GHz, but is consistent with those by Massardi et 
al. (2011) who found that only a few of the 189 sources in 
their complete sample with S20GHZ > 500 mjy have a|[] > 0. 

On the other hand, all SS sources have high-frequency 
spectral indices a|o > — 1 so that their SEDs vS{y) peak 
above the frequency range covered by PACO observations. 
Taking into account the WISE photometry we find that the 
vast majority of them have fp ED < 10 5 GHz and are there- 
fore classified as LSPs, according to Abdo et al. (2010). Their 
median peak wavelength is Ap ED ~ 93 /im. Only 6 of our 
blazars are likely HSPs. 
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